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Abstract 

We analyze the particle-like excitations arising in relativistic field theo- 
ries in states different than the vacuum. The basic properties characterizing 
the quasiparticle propagation are studied using two different complemen- 
tary methods. First we introduce a frequency-based approach, wherein the 
quasiparticle properties are deduced from the spectral analysis of the two- 
point propagators. Second, we put forward a real-time approach, wherein 
the quantum state corresponding to the quasiparticle excitation is explicitly 
constructed, and the time-evolution is followed. Both methods lead to the 
same result: the energy and decay rate of the quasiparticles are determined 
by the real and imaginary parts of the retarded self-energy respectively. Both 
approaches are compared, on the one hand, with the standard field-theoretic 
analysis of particles in the vacuum and, on the other hand, with the mean- 
field-based techniques in general backgrounds. 
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1. Introduction 



In this paper we examine the elementary particle-Uke excitations in generic quantum 
states, as seen from the viewpoint of relativistic quantum field theory. We will refer to 
those excitations as quasiparticles. Quasiparticles are one of the most ubiquitous con- 
cepts in physics, appearing in many different contexts with slightly different meanings, 
for instance in Bose-Einstcin condensation [1], quantum liquids [2], superconductivity [3], 
and, more in general, in condensed matter field theory [4,5], thermal field theory [6,7] 
and AT-body quantum mechanics (non-relativistic field theory) [8] . We shall not by any 
means attempt to review the concept of quasiparticle in this paper. Instead, we shall 
limit ourselves to studying the dynamics of particle-like excitations arising in relativistic 
scalar field theories in states different than the vacuum, focusing on those properties 
which can be extracted from the two-point correlation functions. 

In order to analyze the particle or quasiparticle properties from a second-quantized 
perspective the standard procedure is to add an elementary particle-like excitation to the 
system, and to analyze the subsequent dynamical evolution of the field mode correspond- 
ing to the initial particle momentum. With non-relativistic excitations this procedure is 
clear, given that the particle number is preserved along the evolution and the particle 
concept is well-defined for all times. In this context it can be readily obtained that the 
energy and decay rate of the (quasi) particle excitations can be recovered from the real 
and imaginary parts of the poles of the two-point propagators respectively [8] . 

With relativistic excitations the procedure is less direct since the number of relativistic 
particles in a given field state fluctuates. Properly speaking, in presence of interactions 
particles are only well-defined in the asymptotic limit, where the interaction between 
them can be neglected. Single particle excitations in the vacuum constitute a particular 
case, since Lorentz symmetry is enough to fully characterize their properties [9]. It is 
then a textbook result to show that in the vacuum the energy of a single particle particle 
is also given by the location of the pole of the Feynman propagator, similarly to the 
non-relativistic case. For unstable particles the treatment is less rigorous since properly 
speaking no asymptotic states can be associated to them. Nevertheless by analogy with 
the stable case, and appealing to the optical theorem, it can be additionally argued that 
the imaginary part of the pole of the propagator determines the lifetime of the particle. 

In non-vacuum states, where it is not possible to bring up symmetry considerations, 
it is a priori not completely obvious how the relativistic quasiparticles should be treated 
within a second-quantized formalism. One possibility is to simply forget about quasipar- 
ticles and to study the dynamics of the mean field under arbitrary small perturbations, 
using the linear response theory [6-8, 10]. The response of the mean field is characterized 
by a frequency and a decay rate, which are respectively connected to the real and imag- 
inary parts of the retarded propagator. This constitutes the standard way of analyzing 
relativistic excitations in non- vacuum states [6,8,11], although for typical particle- like 
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excitations the expectation value of the field vanishes, and therefore the excitations con- 
sidered by the linear response method do not correspond to elementary quasiparticles. 
In any case, the approach based on the linear response theory is appropriate to analyze 
another different regime, the hydrodynamic or fluid regime [12-14]. 

Whenever single quasiparticlc excitations are important, a second possibility is try- 
ing to develop an appropriate quasiparticlc framework, following what is done either in 
non-relativistic field theory or in relativistic field theory in the vacuum. This will be 
the main object of analysis of this paper. We will investigate several questions: first, 
which is the field state corresponding to quasiparticlc excitations; second, whether quasi- 
particlc excitations appear naturally in physical situations; third, how to extract their 
properties from the field-theoretic correlation functions, and, finally, what is the range 
of applicability of the quasiparticlc formalism. \^^o will however not discuss under which 
conditions a given field theory in a given field state possesses quasiparticlc excitations: it 
will be simply assumed that the field interaction and state are such that a quasiparticlc 
regime exists. 

There have been some partial results in this direction. Weldon [15] recognized that in 
thermal field theory the imaginary part of the pole of the retarded propagator is linked to 
the decay rate of the elementary excitations. Donoghue et al. [16] discussed the role of the 
the real part of the pole of the propagator as related to the energy of the quasiparticles, 
and the fact that it is not necessarily Lorentz-invariant. Narnhofer et al. [17] studied the 
properties of the quasiparticles in thermal states. Weldon [18] and Chu and Umezawa 
[19] attempted a reformulation of thermal field theories in terms of what they called 
stable quasiparticles, the latter using a thermo-field dynamics formalism. Also in the 
thermo-field dynamics context, Nakawaki et al. [20] studied quasiparticlc collisions (a 
subject which will not be discussed in this paper). Finally, let us mention that Greiner 
and Leupold [21], among other things, studied the conditions for the existence of a 
quasiparticlc regime and analyzed the fluctuations of the quasiparticlc number. 

The kinetic approach to field theory by Calzetta and Hu [22] also goes beyond the 
mean-field treatment. Calzetta and Hu focus on the dynamics of the propagator by 
studying the two-particle irreducible effective action, which leads to equations of motion 
for the correlation functions, from where the equations of motion for the quasiparticlc 
occupation numbers are extracted, thereby connecting the second-quantized formalism 
with a first-quantized statistical description. In this paper we will strictly limit to a 
second-quantized approach; see Ref. [23] for some comments on the connection with the 
first-quantized analysis. 

It should be noted that, since we shall work within a second-quantized formalism, we 
will not study the quasiparticlc themselves but rather the field modes corresponding to 
the initial momentum of the quasiparticles. We will therefore not attempt to follow the 
quasiparticlc trajectory. It will be important to bear in mind that the appearance of a 
dissipative term does not necessarily mean a quasiparticle decay, but most of the times 
it simply indicates a change of momentum of the quasiparticle. 

The paper is organized as follows. In the remaining of this section we introduce some 
of the tools and concepts that will be used in the paper. Of particular interest is the 
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analysis of the two-point propagators in general backgrounds. In Sec.|2]we recall how the 
particle interpretation can be recovered in field theory in the vacuum, introducing the two 
techniques that later on will be employed for the analysis of the quasiparticles. In Sec. [3] 
we develop a spectral approach to the analysis of the quasiparticle excitations in field 
theory, in parallel to the procedure in the vacuum. In Sec.|4]we present a complementary 
real-time analysis of the same problem, namely we study the time evolution of the 
relevant observables in the presence of quasiparticle excitations. In the process we discuss 
the form of the quantum states associated to quasiparticles, and analyze the appearance 
of these states in physical situations. In Sec. [5] we recall the standard mean-field-based 
techniques and compare them with our methods and results. Finally, in Sec. [6] we 
summarize and discuss the main points of the paper. Appendices contain background 
reference material as well as technical details of some calculations. 

Throughout the paper we use a metric with signature (— , -|-, -|-, -|-), work with a system 
of natural units with h = c = 1, denote quantum mechanical operators with a hat, and 
use a volume-dependent normalization in the definition of the field modes [see Eq. (1.17) 
below]. The same symbol will be later used for a quantity and its Fourier transform 
whenever there is no danger of confusion. 



1.1. Quasiparticles in relativistic field theory 

For the purposes of this paper we consider that a quasiparticle is a particle- like excitation 
which travels in some background and which is characterized by the following properties: 

1. It has some characteristic initial energy E. The fluctuations of the energy are 
much smaller than this characteristic value: AE <^\E\. 

2. It has some characteristic initial momentum p. The fluctuations of the momentum 
are much smaller than the characteristic energy: |Ap| <^ \E\. 

3. It has approximately constant energy and momentum during a long period of time 
T = 1/r, before it starts to decay. Here "long" means that the decay rate F has 
to be much smaller than the de Broglie frequency of the quasiparticle: F <C \E\. 

4. It is elementary, meaning that it cannot be decomposed in the (coherent or inco- 
herent) superposition of two or more entities, having each one separately the same 
three properties above. 

Notice that the third property is somewhat redundant, since by the time-energy uncer- 
tainty principle the energy fluctuations are at least given by the decay rate: AE > F. 
The quasiparticle is essentially characterized by the energy E, the momentum p and 
the decay rate F. Besides that, the quasiparticle can be characterized by other quan- 
tum numbers such as the spin (although in this paper we shall only deal with scalar 
quasiparticles). 

If the initial background state has large momentum or energy fluctuations the per- 
turbed state inherits them, and therefore the requirement that the fluctuations of the 
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momentum and energy of the quasiparticles are small cannot be fulfilled. For thermal 
and, more generally, for Gaussian states, we will see that momentum fluctuations are 
comparable to the average momentum when the occupation numbers are of order one. 
Therefore for bosonic systems the quasiparticle description of Gaussian states requires 
relatively small occupation numbers. When the occupation numbers are of order one 
or larger, a quasiparticle description might not be very adequate and it might be more 
useful to move to a hydrodynamic description [12-14] (see Sec. [6] for a further analysis 
of this point) . 

In any case, quasiparticle excitations may exist even in strongly correlated systems. 
The excitations in these systems are usually radically different to the original particles 
which constitute the medium. For instance, the spin and mass of the quasiparticles might 
have nothing to do with the original vacuum particles. The quasiparticles which bear 
no resemblance with the original particles constituting the medium, and which usually 
involve the system as a whole, are also called collective modes [8] . Sometimes in the same 
system the collective modes can coexist with the quasiparticles directly corresponding 
to the medium constituents. 

In this paper we will consider the quasiparticle excitations of a single scalar field (p in 
a given background state. This field can have two different interpretations depending on 
the context. First, the field (j) can have a straightforward interpretation as the fundamen- 
tal field whose excitations in the vacuum give rise to the particles, and whose excitations 
in states different than the vacuum constitute the corresponding quasiparticles. For 
instance, this would be an appropriate interpretation when studying the behavior of a 
pion in a thermal bath. But the field (p can have a second different interpretation as an 
effective field whose excitations correspond to the quasiparticles propagating in the sys- 
tem, not necessarily having any direct correspondence with the fundamental fields. This 
second interpretation is more general and can accommodate for the collective modes. 
For instance, the field (/> could represent the field of sound wave excitations in a Bose- 
Einstein condensate. In any case, in this paper we will consider the scalar field cp as given 
and will not investigate its relation with the fundamental constituents of the system. 

1.2. Dispersion relations 

The dispersion relation is the expression of the energy of the quasiparticle as a function 
of the momentum, namely 



where E is the energy of the quasiparticle excitation, and p is the momentum, with 
the assumption of small spreads. The effective mass is the value of the energy at zero 
momentum, rngfj = Sp=o- When the states are thermal, the effective mass is called the 
thermal mass. For a unstable system with energy E and decay rate F, we define the 
complex generalized energy £ as S'^ := E"^ — iET. The generalized dispersion relation is 
the function giving the generalized energy of a quasiparticle in terms of the momentum: 



E = E, 



P' 



(1.1a) 



£^ = El- iEpTp. 



(1.1b) 
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Notice that the imaginary part of the generahzed dispersion relation places a lower bound 
on the uncertainty of the real part. 

In flat spacetime in the vacuum, the propagation of a stable particle is fully charac- 
terized by the physical mass mp^: 

E^ = ml^ + p'^. (1.2a) 

If the particle is unstable, the generalized dispersion relation is determined by the phys- 
ical mass and the decay rate in the particle rest frame 7: 

"^^ = "^ph + - «"^ph7- (l-2b) 

The decay rate must be much smaller than the particle mass; otherwise one would speak 
of resonances rather than unstable particles. In any case, the generalized dispersion rela- 
tion can be extracted from the location of the poles in the momentum representation of 
the Feynman propagator. In other words, the Feynman propagator can be approximated 
around the particle energy as: 

—iZ 

Gf{uj, p) ^ — 2 ■ ' id'^ ^p"^ + ra\. (1.3) 

-w^ -I- -I- my^ - imp\n ^ 

Several differences with respect to the vacuum case arise when studying the propa- 
gation in a non- vacuum state. First, in general the dispersion relation needs not be 



of any of the forms (1.2), and therefore the effective mass does not determine by itself 
the dispersion relation. The background constitutes a preferred reference frame, thus 
breaking the global Lorentz invariance, and, because of the interaction with the envi- 
ronment modes, the quasiparticle energy is affected by the background in a way such 
that the dispersion relation is no longer Lorentz invariant [24-27]. Moreover, even if 
some of the basic features of the particle propagation can be encoded in the form of a 
dispersion relation, one should bear in mind that a much richer phenomenology appears 
in general (including, for instance, scattering, diffusion and decoherence). Additionally, 
let us point out that it is not completely obvious how in general the dispersion relations 



should be extracted from the poles of a propagator, in a similar way to Eq. (1.3); we 
will also address this point in this paper. 

1.3. Propagators and self-energies in an arbitrary background state 

In the vacuum the analysis of the Feynman propagator is usually sufficient. In a generic 
state p this is not usually the case, and an analysis of the different propagators is in 
order [28-30]. The Feynman propagator, positive and negative Whightman functions 
and Dyson propagator. 



Gii(x,x') = Gy{x^x') 
G\2{x,x) = Gj^{x,x') 
G2i{x,x) = G^{x,x') 
G22{x,x) = Gy){x,x') 



Tr (pr0(x)(^(x')), (1.4a) 

■ Ti{pij){x)^{x')), (1.4b) 

■-Tt{p^{x')^{x)), (1.4c) 

: Tr (pf0(x)0(x)), (1.4d) 
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appear in the closed time path (CTP) formalism (which is natural when dealing with 
non- vacuum states; see appendix |A|, and can be conveniently organized in a 2 x 2 matrix 
Gab, the so-called direct basis: 

^ab[X, X ) - ^^^^^^ ^^^^^ J . (1.5j 

We may also consider the Pauli-Jordan or commutator propagator, 

G{x,x') :=Tt{p[^{x),^{x')]), (1.6a) 

and the Hadamard or anticonmutator function 

G^'^\x,x') := Tt (p{0(x),(^(x')}). (1.6b) 

For linear system^ the Pauli-Jordan propagator is independent of the state and carries 
information about the system dynamics. Finally, one can also consider the retarded and 
advanced propagators, 

Gr(x,x') = e(x°-x'°)G(x,x') = e(x°-x'°)Tr(y5[0(x),0(2;')]), (1.7a) 
Ga{x,x') = e{x'° - x'^)G{x,x') = e{x'^ - x^)Tt {p[^{x),4>{x')]), (1.7b) 

which also do not depend on the the state for linear systems. 

By doing a unitary transformation of the propagator matrix it is possible to work in 
the so-called physical or Keldysh basis [30,31], 

r' ^gW(x,x') Gn{x,x')\ 

Ga'b'{x,x)-^^^^^^^,^ J' ^^■^> 

For linear systems, the off-diagonal components of the Keldysh basis carry information 
on the system dynamics, whereas the non- vanishing diagonal component has information 
on the state of the field. 

The correlation functions in momentum space are defined as the Fourier transform 
of the spacetime correlators with respect to the difference variable A = x — x' keeping 
constant the central point X = (x + x')/2: 

Gabip; X) = j d^x e-"P-^ Gab{X + A/2, X - A/2). (1.9) 

For homogeneous and static backgrounds the Fourier-transformed propagator does not 
depend on the central point X. The retarded and advanced propagators, which are 
purely imaginary in the spacetime representation, develop a real part in the momentum 
representation. 

Obviously not all propagators are independent: the complete set of propagators is de- 
termined by a symmetric and an antisymmetric function. From the Feynman propagator 



""^By linear systems we mean systems whose Heisenberg equations of motion are linear. These correspond 
either to non-interacting systems or to linearly coupled systems. 
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the other Green functions can be derived, but, in contrast, the retarded propagator lacks 
the information about the symmetric part of the correlation function. 

As in the vacuum, self-energies can be introduced for interacting systems. The self- 
energy has a matrix structure and is implicitly defined through the Schwinger-Dyson 
equation: 

Gab{x,x') = + J dSdy G^^X^,y)H^''^y,y')]Gdb{y',z), (i.io) 

where G^^^{x,x') are the propagators of the free theory, and Gab{x,x') are the propa- 
gators of the full interacting theory. The CTP indices a,b,c. . . are either 1 or 2, and 
we use a Einstein summation convention for repeated CTP indices. The ab component 
of the self-energy can be computed in the direct basis, similarly to the vacuum case, as 
the sum of all one-particle irreducible diagrams with amputated external legs that begin 
and end with type a and type b vertices, respectively (cf. appendix [A]) . 

A particularly useful combination is the retarded self-energy, defined as I1r(x, x') = 
T,^^{x,x') + E^^(x,a;'). It is related to the retarded propagator through 

Gn{x,x') = G^^\x,x') + j d^ydSG^l\x,y)[-iT.^{y,y')]GK{y\z), (1.11) 

This equation can be regarded as a consequence of the causality properties of the retarded 
propagator. A similar relation holds between the advanced propagator Gp^{x, x') and 
the advanced self-energy Tix{x, x') = T}^{x, x') + Y?'^{x, x'). Another useful combination 
is the Hadamard self-energy, which is defined as = Y}^[x^x') + 

[or equivalently as Ti'^^\x,x') = —T}'^{x,x') — T?^{x,x')] and which is related to the 
Hadamard propagator through [23] 

= -ij d4ydVGR(x,y)sW(y,y')G'A(y',x') (1.12) 

if the system has a sufficiently dissipative dynamics and interaction is assumed to be 
switched on in the remote past (otherwise the right hand side of the above equation 
would incorporate an extra term) . All self-energy combinations can be determined from 
the knowledge of the Hadamard self-energy and the imaginary part of the retarded self- 
energy. 

So far, all expressions to arbitrary background states p. For static, homogeneous and 
isotropic backgrounds, Eq. ( |1.11[ ) can be solved for the retarded propagator by going to 
the momentum representation: 

— ^ 

G'r(w,p) = 2 I f2 I V / ^ 0^- (^-^^^ 

Notice that in general the self-energy is a separate function of the energy oo and the 3- 
momentum p, and not only a function of the scalar p^, as in the vacuum. The Hadamard 
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function admits the following expression [which can be derived from Eq. (1.12)]: 

iSW(^,p) 



G«(u;,p)=i|GR(u;,p)|2s«(u;,p) 



[-cj2 + El+Re SR(a;, p)]2 + [Im Sr(u;, p)]2 ' 

(1.14) 

From the retarded propagator and the Hadamard function one can show: 

Gf(w,p) = K -, (1.15a) 

[-io^ + El + Re Sr((^, p)] + [Im Sr(^, p)]^ 



, 4-a;2 + E2^ReSR(u;,p)] +iSW(.;,p)/2 

Gd(u;,p) = -, (1.15b) 

[-0.2 + El + Re ER(a;, p)] + [Im SR(a;, p)^^ 



f ^ iSW(c^,p)/2+ImSR(u;,p) 

G_(w,p) = ^ (1.15c) 

[-0.2 + £;2 + Re SR(a;, p)] + [Im SR(a;, p)]^ 

^ iSW(u;,p)/2-ImSR(a;,p) 

G+(u;,p) = K (1.15d) 

[_^2 + + Re 2r(u;, p)] V [Im Sr(.;, p)]=^ 

The imaginary part of the self-energy can be interpreted in terms of the net decay 
rate rp(a;) for an excitation of energy uj — i.e., decay rate T~{uj) minus creation rate 
r+(u;) [15,32]: 

ImER(cu; p) = -cu[r-(c^) - T+{u;)] = -u;Tp{u;) (1.16a) 

In turn the Hadamard self-energy can be interpreted in similar terms as [32]: 

S(i)(u;) = -2i|u;|[rp(Lj) -hr+(Lj)] = -2i\uj\Tp{uj)[l + 2n{uj)], (1.16b) 

where n{uj) is the occupation number of the modes with energy uj. When the field state 
is not exactly homogeneous, the expressions in this paragraph are still correct up to 
order Lp, where L is the relevant inhomogeneity time or length scale. 

In most sections of this paper we discuss general properties of the propagators and 
self-energies, which do not depend on any perturbative expansion. In practice, however, 
many times the only way to evaluate the interacting propagators is through a perturba- 
tive expansion in the coupling constant. For physically reasonable states far ultraviolet 
modes of the field are not occupied. Therefore, there is an energy scale beyond which 
the field can be treated as if it were in the vacuum. In the case of thermal field theory 
this scale is given by the temperature, and the Bose-Einstein function can be viewed as 
a natural soft cutoff to the thermal contributions. The counterterms which renormalize 
the vacuum theory also make the theory ultraviolet-finite in any physically reasonable 
state. This does not mean however that the renormalization process is not modified. As 
we shall see in the following sections, and as it is shown in Refs. [16,32,33], finite parts 
of the counterterms need be adjusted in a different way for each background if a physical 
meaning is to be attributed to the different terms appearing in the action. 
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The infrared behavior is more subtle. On the one hand, naive perturbation theory 
may break down under certain circumstances. In the case of thermal field theory, if the 
relevant contribution to a given process comes from external soft momenta, in order to 
do perturbation theory in a meaningful way the tree level propagators must be replaced 
by Braaten and Pisarski's resummed propagator [34-36], which incorporates the effect 
of the hard thermal loops [6, 37]. On the other hand, even if resummed propagators are 
used, infrared divergences may still arise when the masses of the fields are negligible. 
These divergences may show up in the final results of the calculations, or can be hidden 
in the intermediate stages, leading to finite but incorrect results if they are not properly 
regularized [38]. Additionally, as we shall comment later on, the infrared behavior may 
lead to a modifications to the ordinary quasiparticle decay law [39-42] . The investigation 
of the infrared divergences at finite temperature is still an open problem [37] . 

1.4. Open system viewpoint for the quantum field modes 

The mode corresponding to the propagating quasiparticle can naturally be regarded as 
an open quantum system [43-45]: naively, the field mode would constitute the reduced 
system and all the other modes would form the environment. However, since the mode- 
decomposed field operator is a complex quantity obeying the contraint (f)p = 0!_p, instead 
of focusing on a single mode, it proves more useful to choose as the system of interest 
any two modes with given opposite momentum, and as the environment the remaining 
modes of the field, as well as the modes of any other field in interaction. 
The field (f) can be decomposed in modes according to 

</,p=-^ |d3xe-P-<^(x), (1.17) 

where V is the volume of the space, a formally infinite quantity which plays no role at 
the end. Given a particular momentum p, the system is composed by the two modes 
(j)p and (f>-p, and the environment is composed by the other modes of the field, <^q, with 
q 7^ ibp. Should there be other fields in interaction of any arbitrary spin, the modes of 
these additional fields would also form part of the environment. The entire system is in a 
state p; the state of the reduced system is Ps = Trenv P, and the state of the environment 
is pe = Trsys p. Generally speaking, the state for the entire system is not a factorized 
product state (i.e., p ^ ps® Pc)- 

The action can be decomposed as 5 = Ssys + •S'count + 'S'env + >S'int) where Sgys is the 
renormalized system action. 




(1.18a) 



'S'count is the appropriate counterterm action, 

5count = / dt {(Zp - l)</.p</.-p - [Zp(p2 + ml) - El] (t>p(j>-p] , (1.18b) 
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and S'cnv and •S'int are respectively the environment and interaction actions. Notice that 

1/2 

we have allowed for an arbitrary rescaling of the field (j) (p/^p and for an arbi- 
trary frequency of the two-mode system Ep, which needs not be necessarily of the form 
(p^+m^)^/^. Since it is always possible to freely move finite terms from the system to the 
countertcrm action and vice versa, both the field rescaling and the frequency renormal- 
ization should be taken into account even if no infinities appeared in the perturbative 
calculations. In this paper we will investigate a physical criterion in order to fix the 
values of these two parameters. Notice that this is not just a matter of notation: in 
particular, the election of Ep will determine the value of the energy of the quasiparticle. 

The environment and interaction actions depend on the particular field theory model, 
and in general not many things can be said about them. However, provided that the state 
of the field is stationary, homogeneous and isotropic, under a Gaussian approximation 
the real environment can be always equivalently replaced by a one-dimensional massless 
field and the real interaction can be replaced by an effective linear interaction with 
this environment [32]. In other words, any scalar two-mode system can be equivalently 
represented in terms of a pair of quantum Brownian particles [45-47] , this is to say, by 
a pair of quantum oscillator interacting linearly with a one-dimensional massless field. 
In this paper we will apply the Gaussian approximation, and it will prove useful for us 
to reason in terms of the effective coupling constant and the effective environment. The 
explicit details of the correspondence, which can be found in Ref. [32], will not be needed 
though. 
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We begin by reviewing some aspects of the notion of particle in standard quantum field 
theory in the Minkowski vacuum. While most results in this section can be found in 
standard quantum field theory textbooks (see for instance Refs. [9,48-52]), we present 
them in sonic detail because, first, analogous steps will be followed when studying quasi- 
particle excitations in general backgrounds, and, second, in order to clarify some aspects 
which will prove relevant later on. Except where more detailed references are given, we 
address the reader to the aforementioned textbooks for the remaining of this section. 

Let us consider a scalar field theory whose degree of freedom is the scalar field operator 
(j). If the theory is free, the field operators are connected to the creation and annihilation 
operators via: 

^P = ^=(«P + «-p)> (2-1) 

where = m? + and where a\, and Op are the creation and annihilation operators 
respectively, which verify the following conmutation relations: 

[ap,aq] = [at„ai] = 0, [a^,a\] = = ^{2nf5^''\p - <i), (2.2) 

where we recall that V is the (formally infinite) space volume. If the theory is interacting, 
the same above relations hold in the interaction picture. The normalization is chosen so 
that it closely resembles the quantum mechanical normalization with a finite number of 

degrees of freedom. 

The Hilbert space of the states of the theory has the structure of a Fock space. For 
non-interacing theories, the Fock space can be built with the aid of the creation and anni- 
hilation operators: |pi • • • Pn) = Spi • • • ap„|0) (this equation assumes that all momenta 
are different; if this is not the case, the right hand side should incorporate a factor 1/m! 
for each repeated momentum) . The states this way created have well-defined momentum 
and energy: 

pIpi • • -Pn) = (Pi H h Pn)|Pl • • -Pn), (2.3a) 

^|pi •••?„) = (^(°) + E^, + --- + EpJIpl • • • p„), (2.3b) 

where the momentum and energy operators are: 



kaj^Ok, H = / ^^k aL«k + 7^ • (2.4) 



Any state of the field can be expanded in the Fock space: 

|*) = X] n /(Pl'---'P"i)|Pl---Pm) 



m pi,...,Pm 
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. Complementary to the Fock space expansion, the Hilbert space of the field also admits 
a mode decomposition. Namely, every pure state of the field can be decomposed in the 
following way: |^') = Xlnk Hk '^("'k)|"'k)- where the modes |nk) verify 

p|nk) = nkk|nk), -^|nk) = nk^^kl^k) (2.5) 

The situation gets more involved when the theory becomes interacting. In this case, 
the Hilbert space is still spanned by a set of eigenvectors of the momentum operator 
p and the the full Hamiltonian H. However, multiparticle states cannot be labeled by 
the momentum of each particle in the state, because particles are interacting and the 
momentum of the particles (and even the number of them) changes. However, in the 
remote future and past particles are well separated and the interaction between them 
is negligible. Labeling by |pi • • •Pn)in{out) t^is eigenstate of the full Hamiltonian that 
corresponds to a multiparticle state in the limit t — oo (t — > oo), one has 

-f^IPl • • • Pn)in(out) = (-E'pi H \- -E^p„)|Pl • • • Pn)in(out)> (2-6) 

where in this case = rri^^ + p^, with mph being the physical mass of the particles, 
which differs in general from the bare mass present in the Lagrangian. Notice that the 
particles appearing in the in or out states do not necessarily correspond to the particles 
appearing in the corresponding free theory: any unstable particle will not appear in 
the asymptotic states, and we will possibly have to add bound states to the asymptotic 
states. (For simplicity our notation assumes just one particle species and does not take 
into account these possibilities.) Notice also that the in and out states are defined for all 
times (although they only have special properties in the asymptotic limits). Therefore, 
either the in or out Fock spaces built from those states can be chosen as a basis for the 
Hilbert space of the interacting theory. In the remote past and future one can build a 
free theory that matches the properties of the interacting theory in these regimes. These 
"free states" correspond unitarily to the in and out states of the interacting theory in 
the asymptotic limit^ (the correspondence being different in each case) . The in and out 
states are therefore also unitarily related (via the S matrix). See Refs. [9,51,53] for more 
details. 

2.1. Particles: the spectral approach 

2.1.1. The Kallen-Lehmann spectral representation 

An arbitrary correlation function G^nyix, x') admits the following Kallen-Lehmann spec- 
tral representation, 

/ ds p{s)G^°^ {x,x';s), (2.7) 
Jo 

^In order to properly define this correspondence one must work witii wavepackets, so tiiat tliere is 
localization in time and space; otherwise the states of the Fock space are completely delocalized. See 
Ref. [9] for more details. 
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where G^y{x, x'; s) is the corresponding free function with squared mass parameter s 
and p{s) is the vacuum spectral function, which is defined as 



;lj;(2^)V^)(p-p„)KO|0(O)|a)|^ 
zvr — ^ 



(2. 



and which verifies, among others, the following properties: (i) p(s) > 0, (ii) p{s) = for 
s < and (iii) Jq°° ds p{s) = 1. In this last equation \a) is a complete set of orthonormal 
states, eigenstates of the four-momentum operator p = (If, p), spanning the identity: 



1 = EJ«)(«I, p' 



\a) 



p^\a). 



An example of such class of states are the in or out 



states just considered. 

The spectral representation can be also expressed in momentum space: 

G+(p) = 2^p(V)e(/), 

= 2^p(V)sign(/), 

Jo 

Gf{p) -- 



(p0 + ie)2 + p2 + s' 
-ip{s) ds 



+ s 



(2.9a) 
(2.9b) 

(2.9c) 
(2.9d) 



The first two equations show that in vacuum the Whightman functions and the Pauli- 
Jordan propagator essentially amount to the spectral function. The latter two equations 
show that the retarded propagator and the Feynman propagator have well-defined ana- 
lyticity properties when considered functions in the complex p^ plane. In fact, they also 
have analyticity properties in the complex energy plane, as shown in the following two 
equivalent representations: 



Gr{p) 
Gf{p) 



dk' 



dk^ 



ip{k^^ 



P 



ip{k^^ 



ip{k 



02 



P 



A;0 + ie p^ + k° 



(2.10a) 
(2.10b) 



Taking into account the relation of the spectral function with the Pauli-Jordan propa- 



gator, given by Eq. (2.9b), one can also write 

l-dk^ iG(fcO,p) 

This last equation also follows directly from the definition of the retarded propagator. 



2.1.2. Stable particles 

From all the states of the theory, let us single out the one-particle states corresponding 
to stable particles (assuming they exist), characterized by the momentum p and the 
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physical mass mph: p|p) = p|p), ^|p) = Ep\p), with ii^p = p^ + m'^^. Considering 
those states exphcitly, the spectral function can be developed as 



p{-p') = Z6{p' + m^h) + 0{P' + mt)a{-p'), (2.12) 

where Z = |(0|(/)(0)|p)p is a positive constant (which does not depend on the momentum 
because of the Lorentz invariance of the theory), a{—p'^) is the contribution from the 
multiparticle states, and m^, is the minimum rest energy of the multiparticle states. The 
Z constant is frequently renormalized to one by rescaling the field (which amounts to 
adding a suitable counterterm to the original action). 



According to Eq. (2.9d), the Feynman propagator in momentum space has the struc- 
ture: 

-iZ , -ia(s) 

_|_ _ p"^ + s — le 



G^^P) = + . d.-^^^. (2.13) 



Thus, the Feynman propagator is analytic in the complex p plane, except for a pole 
at = — ?7iph + and a branch cut starting at p'^ = —ml + ie and running parallel 
to the real axis. There are additional poles if bound states can be formed; we shall not 
take into account this possibility in the following. In presence of massless particles the 
pole and the cut may overlap. When doing perturbation theory this is reflected in the 
appearance of infrared divergences in the calculation of the residues of the pole. For 
a discussion of this point in the case of QED see e.g. Ref. [52]. We will assume in the 
following that the pole and the cut are well-separated. 

In perturbative field theory, the renormalized Feynman propagator can be resummed 
in the following way 

pZ _|_ j^Z _|_ ^(^—p^j 

where m is the renormalized mass and S(— p^) is the self-energy., corresponding to the 
sum of all one-particle irreducible Feynman diagrams. The locus of the poles and cuts is 
given by the solution of the equation p^ + m? -|- S(— p^) = 0. The zeros of this equation lie 



the equation is the pole at — 



next to the real axis as dictated by the spectral representation (2.9d). The lowest zero of 



renormalized mass coincides with the physical mass, m = ?n.pj^, so that S(— mpj^) = 0. 

Although it is usual to consider the Feynman propagator, the analytic structure of 
other propagators can be studied as well. The retarded propagator will be of special 
importance: 

Gk{p) = . + r^'^i^^^o- (2-15) 

pz _|_ j^z _ ^gpu p'^ + s — lep^ 



As shown in figure 2.1 in the complex energy plane (complex plane), the retarded 
propagator is analytic except for the two poles at p^ = =bSp — ie and the two branch 
cuts. 
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Imp" ' 




■ • 


• ■ 


-\/ml + p'' 


-V"»ph + P' 


V'"»ph + p' 


^ml + p2 Rep" 



Figure 2.1.: Analytic structure of the retarded propagator in the vacuum, as seen in 
the complex energy plane. There are two poles corresponding to the stable 
particle, and two branch cuts, whose branching points indicate the minimum 
energy for the multiparticle states. Between the poles and the branching 
points there might be as well other poles corresponding to bound states 
(not shown in the plot). 



2.1.3. Unstable particles 

Recall that unstable particles are not asymptotic states and hence do not correspond 

to eigenstatcs of the Hamiltonian nor they belong to the 1-particle sector of the Fock 
space. Rather, they are combinations of those multiparticle states corresponding to the 
decay products of the unstable particle. In terms of the Kallen-Lehmann representation 
of the propagator, this means that unstable states are represented by a branch cut 
rather than a pole. Hence the spectral function for the multiparticle states is given by 
p{—p^) = 0{p^ + m1)p{—p'^), where m* is in this case the minimum energy to create a 
multiparticle state. 

Unstable particles however develop a pole of the propagator in a second Riemann 
sheet of the complex plane [48,54-59]. The Feynman propagator can be resummed 
as 

= p2 + rn? + E(-p2) = p^ + w? + Re S(-p2) + i im S(-p2) " ^^"^^^ 
We now define the mass of the unstable particle as the lowest order solution of 



- rWpi^ + + Re S(mph) = 



and identify 



7 



ImS(m^i^) 



(2.17) 



(2.18) 



as the decay rate in the particle rest frame, according to the optical theorem. The 
mass as defined above corresponds approximately to the rest energy of the particle (this 
assertion will be checked later on), although it should be noted that the energy of an 
unstable particle fluctuates according to the time-energy uncertainty principle. Thus, 
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the second Riemann sheet of the Feynman propagator in momentum space has a pole 
in the region ImS(p^) < 0, whose real part corresponds to the approximate mass of the 
particle and whose imaginary part corresponds to the decay rate: 

Gf{p) ~ 2 ^ ■ + G{p), (2.19) 

The function G{p) is analytic function in the vicinity of the pole, but it does develop a 
singular behavior when approaching the different particle creation thresholds. 

2.2. Particles: the real-time approach 
2.2.1. Time-evolution of the propagators 

So far we have carried out the analysis in the energy-momentum representation. It will 
prove also illustrative to consider the time-momentum representation of the propagator, 

,0 



GF(M';p) = /^e-A*-OGp(p). 



The aim is to compute the behavior of the propagator for large time lapses. We shall 
consider both the stable and unstable cases simultaneously (with 7 = e if the particle is 
stable). 



From Eq. (2.19), the time behavior of the pole can be easily derived: 



GF{t, t'; p) = . = e-Vp'+^^ph-^'"^'^^!*-*'! ^^^^^ 

O /t^S _1_ m'2 /I'm , ^/ 



2^p^ + m^^ - iruph-f 

In order for the particle concept to be meaningful, the condition 7 <C mph must be 
verified (otherwise one would speak of wide resonances rather than particles). Under 
these conditions, the above expression can be approximated as 

Gp{t, t'; p) = =^ e-^i^pl*-*'l e-rpl*-*'l/2 +G{t, t'; p) (2.20) 
2Ep 

with Ep = p^ -|- 'iT'ph, and where we have defined the decay rate in the laboratory rest 
frame 



Tp:='^ = -^lm^{ml^). (2.21) 



The behavior of remaining piece G{t, t'; p) remains to be determined. In general, it 
depends on the precise value of the spectral function along the branch cut, which in 
turn depends on the multiparticle structure of the theory. However, by appealing to the 
Riemann-Lebesgue theorem, one can show [48] that an order of magnitude estimation 
of G{t, t'; p) is given by 

\G{t,t';p)\r^ (2.22) 
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where M is the scale in which the leading multiparticle threshold starts, and where a is 
some positive coefficient which depends on the particular structure of the theory. 

For unstable particles three different time regimes can be therefore distinguished. For 
short times, of the order of the de Broglie size of the particle m~^^, transient effects dom- 
inate, and the behavior depends on the particular details of the field theory model. For 
large time lapses {\t — t'\ » M^^), the behavior is dominated by the pole, which in this 
case is located off the real axis. The modulus of the propagator decreases exponentially 
with a rate Fp. Transient effects are subdominant in this regime, since they decay in 
a much faster timescale M~^. For extremely long times, however, "transient" effects 
dominate again, since any power low decay dominates over an exponential decay for 
sufficiently long times. Experimentally the breakdown of the exponential law for very 
long times is almost never observable, since power-law terms dominate again after many 
particle lifetimes, when the chances of observing the particle are almost null (as we will 
see in the following). 

2.2.2. Two-point functions and asymptotic fields 

We have just seen that the 2-point correlation functions match to the 2-point correlation 
function of a free field plus an additional multiparticle contribution, which vanishes for 
long times: 

G+{t,t';p) = ZG^^''\t,t';p) + G+{t,t';p) — . ZG':l''\t,t';p) (2.23a) 

t— r — »oo 

or equivalently 

(O|^pC/(i,t')0-p|O) — ^ ^ {Pmt,t')\p) (2.23b) 

where |0) and |p) are the vacuum and 1-particle state of the interacting theory respec- 
tively, and where U{t,t/) is the time evolution operator. Therefore, we can make the 
identification 

IP) = ^^-p|0)- (2-24) 

The symbol = here means equivalence when evaluated in a matrix element in the limit 

of large time lapses. Physically, the field operator excites the one-particle state and the 
multiparticle sector, but multiparticle excitations are off the mass shell and they decay 
quickly. 

The heuristic argument given above can be connected to the fact that the particle 

content of the theory corresponds to that of a free theory in the asymptotic limits. Let 
us consider the asymptotic field operator (p, which by assumption obeys free equations 
of motion, 

id^d'' + ml^)^ = 0, -|p + (p2 + J^p = 0, (2.25a) 
and which corresponds to the field operator through 

0(x) ^ Z-^/2^(x), ^p^Z-^/^p- (2.25b) 
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Let us also consider the corresponding creation and annihilation operators, 



(2.26) 



where VTp is the canonical momentum operator and E 



ph- 



In and out Fock 



spaces can be constructed with the field [9,51,53]|^ With the asymptotic field operator 
the one-particle state can be represented as 



ip) ^ y2^^-p|o) 



in/out 



st,|o) 



in/out 



(2.27) 



Comparing Eq. (2.24) and Eq. (2.27) we see that the distinction between the field </> and 



the the asymptotic field (\) may be blurred provided the constant Z is renormalized to 
one and assumes that in the asymptotic limit the field obeys efi'ective free equations of 
motion. Under these assumptions one may simply write 



ip) ^ y2^0_p|o) 



in / out 



in/out 



(2.28) 



When particles are unstable the situation is less clear since strictly speaking a one- 
particle sector which also is an eigenstate of the Hamiltonian does not exist (the only 
eigenstates being the multiparticle states corresponding to the decay products of the 
unstable particles). However, if particles are long-lived one may think of approximate 
1-particle states (which in fact correspond to multiparticle state combinations). We 



assume that those are also approximately given by Eq. (2.28). 



2.2.3. Physical interpretation of the time representation 

We have previously seen that in frequency space the poles of the propagator are con- 
nected with the energy and lifetime of the particle. Let us confirm this statement study- 
ing the time evolution of the expectation value of the Hamiltonian operator. In field 
theory, strictly speaking, only asymptotic properties are completely well-defined. How- 
ever, as long as sufficiently large time lapses are considered (compared to the relevant 
inverse energy scales), approximate results for finite times can be obtained. We shall use 
the asymptotic particle representation in terms of fields (described in the previous sub- 
section) for finite but sufficiently long time lapses, shall blurry the distinction between 
the usual and asymptotic fields (supposing that the quantity Z has been renormalized 
to one) and shall treat the stable and unstable cases simultaneously. 



According to Eq. (1.18a), the reduced Hamiltonian operator for the relevant two-pair 
mode is given by 



sys 



(2.29) 



^One possible concern is the fact that there is an apparent contradiction between the two canonical 
conmutation relations \c\>, vr] = i and [(^, 7f] = i. In reality there is no such contradiction because 
the two fields are equivalent in the asymptotic limits only when evaluated in matrix elements (weak 
operator convergence) [51]. In any case, the reader must be warned that a fully rigorous analysis is 
usually not feasible for standard field theories. 
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If in the remote past time to a particle is introduced into the system, the state for later 
times will be |p) = ap|0). The evolution of the number of particles in this state is given 
by 



E{t,to;p) := (p|^p(t)|p) = ^p(0|ap[/(to,t) (at Sp + alpa_p + ij [/(t, to)«J,|0). 

(2.30) 

Introducing a resolution of the identity we obtain: 

E{t, to; p) = Ep J2{0\apU{to, t)al,\a){a\apU{t, to)«J,|0) + E^^\ 

a 

where = (0|-^sys|0) is the vacuum energy. By energy and momentum conservation, 
only the vacuum survives from the above summation. Therefore: 

E{t, to; p) = + 4El\{0\4>-pU{to, t)<Ap|0)P, 
from where we obtain 

E{t,to;p) = e'^^') +4El\GF{to,t;p)\^. (2.31) 



Introducing the explicit value of the propagator, given by Eq. (2.20), and neglecting the 
off-shell contribution G{t,t';p) we get the expected result 

E{t, to; p) = + Ep e-rp(*-*o) . (2.32) 

Particles have energy Ep and decay in a timescale Fp in the e domain of validity of the 
exponential law. 

We recall once more that unstable particles do not correspond to any eigenstate of 
the Hamiltonian, and thus no asymptotic states can be associated to them. Therefore, 
strictly speaking, one-particle states are a linear combination of many multiparticle states 
and one cannot make reference to energy conservation, since energy conservation is asso- 
ciated to asymptotic properties. However, if the lifetime of the particles is long enough 
one can think of approximate asymptotic states and approximate energy conservation, 
so that the above calculation would be approximately valid. 
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3.1. Kallen-Lehmann spectral representation 

In a general background the quickest and clearest way to derive the spectral repre- 



sentation is by simply recalling the definition of the propagators. From Eq. (1.7a) we 
immediately obtain the spectral representation for the retarded propagator [28] 

GKiP;X)=J ^^y^y^. (3.1a) 



This equation is identical to its vacuum counterpart, Eq. (2.11). The advanced propa- 
gator follows a similar representation, 

/-dfc" zG(fc°,p;X) 

Ga{p;X) = J— ^^_^^_^^ . (3.1b) 

Hereafter the Pauli-Jordan function will also be called spectral function. The similarities 
with the vacuum case end here. The retarded and advanced propagators are the only 
propagators that can be expressed in terms of an integral of the spectral function. Notice 
also that in general the spectral representation can only be expressed as an integral over 
the energy, and not as an integral over the invariant mass. 
The Pauli-Jordan function verifies the following properties: 

G{p;X)>0, if/>0, (3.2a) 
G{-p;X) = -G{-p;X), (3.2b) 

^A;°G(fc;X) = 1. (3.2c) 

The first two properties are a simple consequence of the definition of the propagator. 
The third property is a sum rule, consequence of the equal-time commutation relations, 
m,^),dt^{t,y)]=i6^^\^-y). 
For stationary backgrounds an explicit representation for the Pauli-Jordan function 



can be obtained, similarly to Eq. (2.8). The Pauli-Jordan function can be expressed in 



the basis of eigenstates of the Hamiltonian as 

G{t,t';p) = ^p„(a|[<^_p(t),^p(0]|a) = ^Pa{a\4>-p e-^(*-*') c^p e*^(*-*') \a) - (c.c), 

a a 

where pa = {a\p\a). We have used the fact that the state is stationary, so that the 
density matrix operators diagonal in the basis of eigenstates of the Hamiltonian. Let us 
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now introduce a resolution of the identity 1 = 

G(i,t';p) = 5]pa(a|0_p e-^^(*-*') |/3)(/?|^p e^^(*-*') |a) - (c.c.) 

a,f3 

= Y,Pa{a\^~pmmp\a) e-(^/5-^")(*-*') -(c.c.) 

a,f3 

= Y,Pa\{a\^pm\-2i) sin [{Ef, - E^){t - t% 

where we have used that = 0_p Going to the frequency space we obtain the desired 
expression: 

G{p) = J2 Pa\{a\^p\P)f [S{p'' + E^- Ep) - 5{p^ -E^ + Ep)] . (3.3) 

For stationary background states the Pauli-Jordan propagator is proportional to the 
probability for the field operator of momentum p to induce a transition to a state with 
higher energy p^^ minus the probability to induce a transition to a state of lower energy 



3.2. Spectral analysis of the quasiparticles 

Let us assume that around some range of energies the Pauli-Jordan propagator has the 
following structure: 

G{u, P) = ^ [^{^ - Ep) - 6{uj + Ep)], for ~ Ep. (3.4) 

This means that the field operator with momentum p creates an excitation whose energy 

is exactly Ep. Since there is no spread in the energy, the excitation must be infinitely 
lived (must be stable). The corresponding retarded propagator is: 

GR{io,p)= "'^P . , for |cc;| ~ £;p. (3.5) 
— a;^ + E^ — lue 

The excitation this way created would have exact energy Ep and exact momentum p. 
Therefore it corresponds to a stable quasiparticle. Notice that, in contrast to the stable 
particles in the vacuum, Zp can depend on the 3-momentum p and E^ need not be of 
the form p^ -|- w?. 

However, stable quasiparticles are an idealization, and do not correspond exactly to 
any physical situation, since dissipation is a generic feature of non-vacuum states as 
long as there is interaction. We recall that by "dissipation" we do not necessarily mean 
the quasiparticle decaying into a product of different quasiparticles: "dissipation" can 
simply mean the quasiparticle changing its momentum. In any case, stable quasiparticles 
may be a good approximation in many situations in which the effective coupling to the 
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environment is very small (in the sense of the quantum Brownian motion correspondence 
explained in Ref. [32] and summarized in the introduction). However, this does not mean 
at all that the real coupling must be small, or that the approximation is limited to weakly 
interacting systems; as a matter of fact, in strongly coupled situations there might be 
situations in which assuming free quasiparticles might well be a good approximation. 

Anyway, quasiparticles are in general not stable, and instead decay with some rate 
Fp. The generic form of the retarded propagator in terms of the self-energy is 

— 1 

Gr{uj,p) = — — . (3.6) 

-uj-^ + + Sr(w,p) 



Following the vacuum analysis of subsection 2.1.3 whenever there is a long-lived quasi- 



particle the analytic continuation of the retarded propagator can be approximated by: 

Gk{uj, p) = 2 ■ r + P)' (3-^) 

where G{ijJ, p) is an analytic function in the vicinity of the pole. This means that for 
unstable (but long-lived) quasiparticles, the Pauli-Jordan function, instead of being a 
delta function around the energy of the pole, it is an approximate Lorentzian function, 
whose width corresponds to the decay rate. The momentum-dependent functions Ep 
and Tp are given, in terms of the self-energy, by 

=m2 + p2+ReSR(£;p,p), (3.8a) 

rp = --^ImSR(ii;p,p), (3.8b) 
Ep 

under the assumption that Tp is much smaller than E^,. The spectral representation 
indicates that Ep corresponds to the quasiparticle energy. The physical meaning of 
Tp can be extracted from the interpretation of the imaginary part of the self-energy. 



Eq. (1.16a): it corresponds to the net decay rate of the quasiparticle. 



The structure of the retarded propagator (3.7) is completely analogous to that of the 



vacuum (see Sec. [2] and figure 3.1). There are however two important differences that it 



is worth commenting. First, in the vacuum the analysis can be equivalently performed 
with the Feynman or the retarded propagators, while in generic backgrounds the spectral 
analysis can only be applied to the retarded propagator, since in general the Feynman 
propagator has an additional explicit dependence on the background state of the field. 
Second, in the vacuum one can study the spectral structure either in terms of the energy 
or in terms of the squared four-momentum. In general, there is a preferred reference 
frame and therefore the explicit Lorentz invariance of the results is broken, and, hence, 
only a spectral analysis based on the energy is meaningful. See Refs. [11,60] for a more in 
depth investigation of the analytic structure of the propagators in thermal field theory. 

It must be noted that in presence of massless particles, there are situations in which 
the decay rate appears to be infrared-divergent from a perturbative calculation. This can 
be attributed to the fact that in presence of massless particles the retarded propagator 
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Figure 3.1.: Analytic structure of the retarded propagator for quasiparticles in general 
backgrounds, in the complex energy representation. The retarded propaga- 
tor has no poles but only has a cut parallel to the real axis. Its analytic 
continuation (or second Riemann sheet) has a pair of poles, whose real part 
corresponds to the energy of the excitation and whose imaginary part cor- 
responds to the decay rate. The cut need not extend from —00 to -|-oo: it 
might be interrupted for some energy sectors. 



does not necessarily exhibit the form (3.7), because there is no threshold for the creation 
of massless excitations. For instance, when computing the lifetime of a quasiparticle in 
a very hot QED plasma it is found that the retarded propagator does not exhibit any 
singularity near the particle resonance energy [40,41,60], although the propagator is still 
strongly peaked around this point. Although we will comment further on this possibility, 
in the following we will limit our calculations to those situations in which a well-defined 
decay rate can be associated to the quasiparticles. 

When the retarded and advanced propagators can be approximated by 



G'R(a;,p) 
Ga(w,p) 



-uj^ + Rl- iwFj 



the Pauli-Jordan propagator (or spectral function) is given by 



(3.9a) 
(3.9b) 



gW( 



UJ,p) 



2\uj\ZpTp 



2 + Eiy + {u;Tp) 



(3.9c) 



-co 



So far we have analyzed the propagators which are independent of the occupation num- 
ber. Other propagators depend on the occupation number. Approximating n(|u;|) by 
the value at the pole, rip := n{Ep), we obtain from Eq. (1.16b) 



S«(£;p,p) « -2ilm^n{Ep,p){l+2np), 



(3.9d) 
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where Up can be interpreted as the occupation number of the mode with momentum p. 
Combining the above equation with Eqs. (1.15) yields 

G(a;,p)p 



2ujZpTp{l + 2np) 



(3.9e) 



Finally, the Feynman and Dyson propagators and the Whightman functions can be 
represented by 



Gf(w,p) 

Gd(w,p) 
G+{uj,p) 



-iZp{-uj^ + E^) + l^lZpFp (2np + 1) 



-to- 



iZpi-LO^ + El) + Ic^lZpFp (2np + 1) 

{-u^2 ^ Elf + {uTpf 
2|a;|Zprp [np + ejiv)] 

+ ElY + {u^Tpf 
2|cj|Zprp [np + 9{-oj)] 
{-u;2 + El)' + {uTpf 



(3.9f) 
(3.9g) 
(3.9h) 
(3.9i) 
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4.1. Quantum states corresponding to quasiparticles: the free case 

Even in absence of interaction it is not trivial to construct a quantum state which verifies 
the quasiparticle properties mentioned in the introduction. As we have seen, in a non- 
interacting theory in the vacuum the one-particle state is naturally represented by the 
action of the creation operator on the vacuum, |p) = ap|0), or equivalently by the action 
of the field operator: |p) = y^2Ep(p^p\0) . Over a homogeneous and stationary state p, 
the positive-energy quasiparticle state can be represented by the action of the creation 
operator: 

'^p' •= ^ \, °P^Qp = ^^11 Pn,{m}{n+l)\{n + l)p,{m}){{n+l)p,{m}\, 



n,{'m\ 



(4.1) 



where the second equality is written in a highly schematic notation in order to avoid 
cumbersome expressions. Here 

np = Tr(pat,ap) (4.2) 

represents the mean occupation number of the mode with momentum p, |np, {m}) is 
a basis of the Hilbert space with the p sector singled out, and Pn,{m} is the diagonal 
representation of the background state in this basis. Notice that the state is normalized: 
TV /5^' = 1. The factors n -|- 1 on the second equality are a consequence of the Bose- 
Einstein statistics: the probability to create an additional particle increases with the 
number of already-existing particles. 

The Bose-Einstein statistics are responsible for the following surprising fact: it is a 
simple exercise to show that the expectation value of the particle number, 

(iVp)(+' = Tr (pWiVp) = Tr (p[,+'at Sp) = -^{{Np + 1)^), (4.3) 

Tip -\- I 

is actually increased more than one with respect to the unperturbed value: 

{NpY^^ =np + l + —^=:np + N^+\ (4.4) 
i + np 

where = 1 -|- 6np/{l + Up) is the number of excitations and 

5nl = {{Np - npf) = {N^) - > (4.5) 

is the dispersion of the number of particles in the background state. For a particle 
number eigenstate dn^ = and A^p'*'' = 1, and for a Gaussian state (see appendix [b| 
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Up + l) and A'^^^' = l + rip. Therefore, when the background is in a state which 



5nl = np 

is not an eigenstate of the Hamiltonian, pp^ represents shghtly more than one additional 
quasiparticle. The reason for this is purely statistical: the highly occupied components 
of the state are more likely to become excited, and therefore they tend to gain statistical 
weight, thereby increasing the particle number. In other words, the action of the creation 
operator has two simultaneous effects: on the one hand, adding a quasiparticle to the 
system; on the other hand, increasing the statistical weight of the highly excited states. 
The statistical contribution is significant when the occupation numbers are large. 

The expectation value of the momentum operator is also affected by this statistical 
effect: 

(P)(+) = TV {pM p) = Tr p at,ap) = piV^^ (4.6a) 
Anyway the momentum per excitation is p. The energy is similarly affected: 



(4ys)^+^ = TV (/5<+^i/sys) = TV 



+ = E^^^ + EqpN^+\ (4.6b) 



where = Tr {pHsys) is the energy of the background and Eqp = Ep = + m? is 
energy per excitation. 

Quasiparticles require small momentum spreads. The spread of the momentum per 
quasiparticle in the case of Gaussian states is 



(Ap)<+' 



QP 



1/2 (4np + 5n^)^/2 

P 1 i ; 
l + Up 



(4.7) 



which is a small quantity if the occupation numbers are small. Likewise, the spreads in 
the particle number and the energy have the same corresponding values. Therefore, in 
order for the state pp ^ to adequately represent a quasiparticle occupation numbers must 
be small as compared to one. 

Since the statistical contribution looks awkward, one can imagine a different way of 
constructing the quasiparticle excitation: 



Pp — Pn,{m}\{n + l)p,{m}){{n+l)p,{m}\. 

n,{m} 



(4.8) 



This state has essentially identical properties to Pp"', except that it contains exactly 
one additional particle. Therefore it looks like that it corresponds more closely to the 
quasiparticle concept we discussed before. However, notice that the above state cannot 
be easily created from the background via the action of the creation and annihilation 



operators. We shall argue in the next section that the state pp , and not pp , appears 
naturally when studying quasiparticle creation processes. 

We shall consider that the state p^' represents exactly one additional real quasiparticle 
of momentum p and energy Ep. The additional contribution to the particle expectation 
number will be called statistical quasiparticle contribution. The statistical quasiparticle 
contribution can be interpreted as being a consequence of the increased knowledge of 
the background state which is derived from the creation of a real quasiparticle. 
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Another novelty is that negative energy excitations, or holes, can also be defined, 
represented by the action of the annihilation operator: 

' ■= ;r"-P^^-P = ;r X] Pn,{m}n\{n-l)-p,{m}){{n-l)^p,{m}\. ug) 

This hole state is also affected by the statistical considerations described above, as it is 
manifest by showing the expectation value of the number operator: 

(ATp) = IV (p(,-)Ar) = np - 1 + ^ =: np - %-\ (4.10) 

Hp 

where A'^^"^ = 1 — dn'p/np is the expected number of negative-energy excitations. There- 
fore the hole state contains at most one negative-energy excitation. This state has 
momentum p, (P)^~^ = pAp, and has an additional amount of negative energy Ep, 

(py-^ = e(^') - EpN^r^ 

For a particle number cigcnstatc Ap"' = 1, but a is that for a Gaussian state the 
number of negative energy excitations is actually negative Ap"' = —Hp. The reason for 
this last surprising fact is that the annihilation operator, besides explicitly removing one 
particle from the state, also enhances the probability of the highly populated sectors of 
the state, and the last effect is actually dominating. We will consider that p'p ^ represents 
a single real hole plus some contribution of statistical quasiparticles of opposite momen- 
tum. Again, the statistical quasiparticle contribution can be thought as a consequence 
of the increased knowledge of the background state coming from the absorption of a 
quasiparticle (or creation of a hole). In any case, the fact that for Gaussian states the 
statistical contribution always dominates is an indication that holes cannot be consid- 
ered true quasiparticles in bosonic systems. A confirmation of this fact is given by the 
spread in the momentum per quasiparticle. 



{^p)'qI ■■= ^ {(P')'-' - [(P)<-f } = |P|(5 + l/r^p)^/^ (4.11) 

which is of the order of the expectation value of the momentum per quasiparticle, even 
if the occupation numbers are small. Even if holes are not true quasiparticles in bosonic 
systems, it will be important to remember that negative energy excitations are possible 
when the states are different than the vacuum. 

So far we have seen that quasiparticles and holes are respectively created by the 
creation and annihilation operators. The field operator (f)p, being a linear combination 
of creation and annihilation operators, creates a coherent superposition of quasiparticles 
and holes. In effect, the state cj)-pp<pp corresponds to the linear superposition of a 
quasiparticle and a hole: 

^-pP^p = ^ («p + «-p) P («P + «-p) • (4-12) 

For states characterized by low occupation numbers the quasiparticle contribution dom- 
inates over the hole contribution. 
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4.2. Creation of quasiparticles 



Let us argue that only the states created with the creation and annihilation operators, 
or equivalently with the field operator, are physically meaningful. In a realistic situation 
quasiparticles are created with the interaction of the field with some external agent. This 
external agent can be modeled by a "quasiparticle emitter" — which essentially coincides 
with the usual particle detector considered in the analysis of the Unruh effect [61,62]: an 
external device described by some quantum mechanical degree of freedom Q (which can 
correspond for instance to a harmonic oscillator or a two- level system). The device starts 
in an excited state. For simplicity, let us assume that the emitter is linearly coupled with 
one pair of field modes, gQQ{t)[(j)p{t) -|-0_p(t)]j^ The coupling constant gg is very small 
so that the external device, other than emitting quasiparticles, does not significantly 
perturb the dynamics of the field. The initial state of the field plus detector is assumed 
to be p |1q)(1q|. The aim is to find the final state for the field p' when the detector 
is measured in its unexcited state |0). 

The time evolution of the entire system under the interaction is given by 



/Ototai(i) = Te 



p<^\lQ){lQ\Te 



where the subindex I indicates interaction picture, and 
coupling is small, the above equation can be expanded as 

ft 



+ 



-p- 



(4.13) 
Since the 



1 - M 9QMs)Qiis)ds 



to 



Ptotal{t) 

or, developing the expression 
Ptotai(i) = p(S>\1q){1q 



|1q)(1c 



1 + i 



{s)Qi{s)ds 



to 



to 



dsgQ[4>i{s),p] ® [Qi(s),|1q)(1q|] 



+ 9q 



to 



(4.14) 



ds / ds'Ms)pHs')(^Qiis)\lQ){lQ\Qi{s'). 



to 



The final state for the system, if the emitter is found unexcited at time t, is given by 
the projection into the ground state of the emitter [59] : 



Tremitter (Ptotal(^) I Oq) (Oq I ) 
TVtotal(/5total(t)|0Q)(0Q|) ' 



(4.15) 



Developing Eq. (4.14) we find 

p' = N ['ds [\s' Ms)p4>i{s') ^ {0\Qi{s)\1q){1q\Qi{s')\0) 



to 



to 



to 



N' ds ds' e'^^'-''^l{s)p4>l{s'), 



to 



^Notice that in this simphfied model the emitter cannot couple to a single mode of the field because 
that couping would not be momentum-conserving (or, from another point of view, would not be 
Hermitian) . 
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where $7 is the frequency of the harmonic oscillator and and A^' are normalization 
constants. Expanding in terms of creation and annihilation operators we get: 



P 



' = N" fds f ds' e^^(^-^') (alp hp g-^^pC^-^') +a\,palp e-^^p(^+«') 

Jtn Jtn ^ 



iEp{s-s') 



+ alp/5 at e-^p('^+-') +appa_p e^^p(^+^') +appat, e^^p(-^') 

Let us assume that the frequency of the oscillator is tuned so that 0, = Ep. In this 
case, for sufficiently large time lapses the dominant value of the integral is given by the 
stationary value of the integral of the first term, which amounts to considering energy 
conservation. In that case 

p' « iV'"(at,pGp + alppa^p) = ^{p^^^ + (4.16) 

Therefore, upon deexcitation of the emitter, the system gets promoted to a superposition 
of two quasiparticle states. The argument could be repeated with the same measuring 
device in the ground state, now interpreted as a particle detector. When the measuring 
device gets excited, the state of the field colapses to a superposition of the hole states 
Pp"' and p'-'p. 

4.3. Quantum states corresponding to quasiparticles: the interacting case 

We have just seen that when no interaction is present, positive-energy quasiparticles are 
created by the action of the creation operator, and negative-energy holes are created 
by the annihilation operator. If interaction is present and quasiparticles are stable. 



similarly to the vacuum case (see section 2.2.2), one can think of defining asymptotic 



quasiparticle fields and states, from which the two-point correlation functions can be 
reproduced. Namely, quasiparticle states correspond to 



1 

and negative energy holes correspond to 



fir, -r r 



/5[r' = ^ a^ppalp, (4.17b) 

Tip 

where np = Tr (papOp) is the expectation value of the number of asymptotic quasipar- 
ticles in the state with momentum p, and we recall the symbol = means equivalence 
when evaluated in matrix elements in the asymptotic limit. The asymptotic creation 
and annihilation operators a and are related to the asymptotic field operator through 



p + alp), (4.18) 
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where Ep is the quasiparticle energy. The asymptotic field, which obeys free equations 
of motion, (p + (m? + Ep)(j) = 0, is connected to the usual field through 

Ip ^ Zp '/'(Ap, (4.19) 



where Zp is defined in Eq. (3.4). See Refs. [17,21] for comparable approaches. 

However, since quasiparticles are generically unstable and therefore not completely 
well-defined from a strict point of view, it is not worth pursuing a very formal description. 
Therefore we will blurry the distinction between the usual and the asymptotic fields, and 
work with the usual field operator but rescaled a factor Zp i.e., we will set Zp = Zp 



see Eqs. (1.18)]. With this assumption equations (4.17) can be more simply restated: 



' - ^rn ^p^^p' ^p'' ^ «-pp«-p- (4-20) 

"-p ~r J- "-p 

Bear in mind that this representation is only valid when studying (approximately) 
asymptotic properties. 

When we presented the open quantum system analysis of the field modes (see the 
introduction), we also introduced two free parameters, £'p and Zp. At this point we 
have, on the one hand, a criterion to fix their value, and, on the other hand, a physical 
interpretation for both of them. With respect to the former, Ep must be fixed to the 
value of the real part of the pole of the propagator, and represents the physical energy 
of the quasiparticle excitation. With respect to the latter, Zp measures the probability 
that the interacting field operator excites the quasiparticle state. Roughly speaking, by 
rescaling the field a factor Zp = Zp we ensure that creation and annihilation operators 
create and destroy quasiparticles with the proper normalization. From now on we shall 
assume that such a scaling has been performed so that Zp = Zp. From a practical point 



of view this amounts to setting Zp = 1 in equations (3.9). 



4.4. Time evolution of the quasiparticle excitations 

The aim is to find the time evolution of the expectation value of the Hamiltonian of 
the system. We shall make use of the open quantum system viewpoint presented in 
the introduction. Let us start with the positive energy excitations, and focus on the 
state with momentum p. The expectation value of the Hamiltonian operator in such a 
quantum system is given by 

E^+\t, to; p) := (i?.ys(t))(+) = TV {alpapU{to,t)H,ysU{t, to)) , (4.21) 

Up + i 



where the system Hamiltonian is given by Eq. ( |2.29 ). Discarding the connected part 



of the 4-point correlation functions, as it is manifest in the following application of the 
Wick theorem [see appendix [B] and in particular Eq. ( |B.8[ )]: 

£■'+'(*, to; P) ~ VtI '^"y" (Psflpflp) Trsys (Ps^sys) 

Tip -\- 1 ^ 

+ EpTT,y, [p,apU{to,t)al,U{t,to)] Tr,y, [p,al,U{t,to)apU{to,t)]Y 
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which can be rewritten as 



Er. 



np + 1 



Trsys [psapU{tQ,t)ayU{t,to)\ 



where E^^^ := Tr^ys (ps-^sys) = E-p{n-p+n^-p+l /2) is the energy of the reduced subsystem 
before the perturbation, or as 



np + 1 



2E, 



-{Ep + idt){Ep - idt,)G+{t, to; p) 



(4.22) 



where we have used the expression of the creation and annihilation operators in terms 
of the asymptotic field and asymptotic canonical momentum, Eq. (2.26), and the "field" 

representation of the canonical momentum, TTp^t) = (j)-p{t). 



Using translation invariance the energy of the perturbation can be rewritten as 



E^+\tM,v)^E^^^ + ^^\I{t,to;vt 
np + 1 



where 



I{t,tQ]Y>) 



1 



2Er 



iEp + idtfG+{t,to;p) 



Using Eq. (3.9h), the expression I{t,to;p) can be written as: 

du 



\uj\iu; + EpyTp [np + e{u;)] 



2tt 



-u>' 



+ El,y + {u;rp) 



(4.23a) 
(4.23b) 

(4.23c) 



The above integral is evaluated by integration in the complex plane in appendix [C| 
Neglecting Fp in front of Ep it is found I{t, to; p) ^ Ep{np + 1) e~^p(*~*o), and therefore 
the result is 

E^+\t, to; p) ~ E^^^ + £;p(np + 1) e-rp(*-*o) . (4.24) 

The factor np + 1 is the statistical factor (recall that for a Gaussian state the expected 
number of excitations is A^p^' = np + 1). Discounting this statistical factor the energy 
of the quasiparticle thus evolves as: 



i?gj^(t,to;p)«i?pe-r-(*~*«). 

For the hole excitations we find a similar result 

E^-\t, to; p) ~ + ^pnp e-rp(*-*o) 



(4.25) 



(4.26) 



Recall that, because of statistical effects, for Gaussian states the initial expectation value 
for the energy of the hole excitations is given by EpUp. 

Equations (4.24) and (4.25) are the main results of this section. From these equations 
Ep can be identified as the energy of the quasiparticles, in agreement to the results of 
the previous section, and Fp as the decay rate of the quasiparticles, according to the 
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interpretation of the imaginary part of the self-energy as the decay rate — recall that 
and EpTp correspond to the real and imaginary parts of the self-energy respectively. 

Notice that four basic assumptions have been used in order to reach these results. 
First, we have introduced an approximately asymptotic representation for the field op- 
erators in terms of the creation and annihilation operators, which is valid if we are con- 
sidering the evolution over times much longer than the characteristic interaction time. 
Second, and related to this, we have used a near-on shell approximation for the propaga- 
tors, which is also correct if we are investigating long times (but not extremely long, as 
discussed in section 2.2). Third, we introduced a Gaussian approximation, which makes 
the problem solvable in term of the two-point correlation functions. Finally, we have also 
assumed that the retarded propagator has the analytic structure given by Eq. (3.7). As 
we have already mentioned, in presence of massless excitations there can be situations 
in which the analytic structure is different, leading to non-exponential modifications of 
the decay law (4.24), as in the case of hot QED plasmas [40-42,63,64]. 

It could be argued that the derivation presented in this section is somewhat cyclic, 
because we are assuming from the beginning that Ep is the frequency corresponding to 
the 2-mode system. Following what is done in Ref. [23], one could instead assume that 
the frequency of the system has an unknown value E'^. We would then find a rapid 
oscillatory behavior for the expectation value of the energy. Only when E'^ = Ep the 
energy follows a smooth exponential decay as is expected on physical grounds. 

We have computed the time evolution of the propagators [and hence the time evolution 
of the energy, according to Eqs. (4.23a) and (4.23b)] by going to the frequency space 
and taking profit of the spectral analysis of the propagators carried out in the previous 
section. Alternatively, it is possible to bypass the spectral analysis by considering the 
equation of motion followed by the propagators, which in general is an integro-differential 
equation. Making a local approximation to the self-energy in terms of the energy an the 
decay rate, the equation of motion turns into an ordinary differential equation which 
can readly be solved [65]. Both methods are physically equivalent and lead to identical 
results. 
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5.1. Linear response theory 

The standard way of analyzing the quasiparticle properties is with the aid of the hnear 
response theory [6-8, 10]. In hnear response theory the Hamiltonian H of a quantum 
mechanical system is perturbed with some external perturbation V{t) at time to, so that 
the total Hamiltonian is H + V{t) for t > to. Given some quantum operator 0{t), it is a 
simple exercise to show that to first order in the potential the expectation value of that 
operator in is given by 

POO 

{dit)) = {dit))o + 6{d{t)) = {dit))o-i / Tr{p9{t-t')[di{t),Viis)]), (5.1) 

J to 

where I indicates interaction picture with respect to the external perturbation and (•)o 
indicates the expectation value in absence of the external potential. 

Linear response theory is usually applied to a free or interacting scalar field theory, 
with the following identifications: 



V{t) = J d^xj(t,x)(/>(x), 0{t) = 0(x). 
Then, from Eq. (5.1), one obtains: 

6{4>{t,^)) = -i [ d''xGn{x,x')j{x'). (5.2) 



For the case of an impulsive perturbation, j{x) = j{q) e^"^'^ 5{t), if the retarded propa- 
gator is approximated as 

iZ 

Gyi{uj,p)^-—— ^ , ujr^Ep, 

{2Ep){uj - Ep + tTp/2) 

the following result for the dynamics of the expectation value of the scalar field is ob- 
tained [6]: 

S{4>it, x)) ^ -'-^e{t) e^(q— ^p*) e-r*/2 . (5.3) 
2Ep 

The fluctuation oscillates with a frequency Ep and decays with a rate rp/2. Those are 
interpreted as the energy and damping rates of the excitations. Notice that the energy 
of the oscillation decays with a rate Tp. 

In non-relativistic A^-body theory, linear response is usually applied to the fermion 
density, instead of the field itself, and the final result depends on the density correlation 
functions [8]. Beyond the field theory context, linear response has many different appli- 
cations in statistical mechanics [10], solid state physics [66], and even gravitation [67]. 
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5.2. Effective dynamics of the mean field 

Another possible method to study the quasiparticle properties has been the analysis of 
the effective dynamics of the mean field. This method has been previously used in the 
literature by Weldon [18] and by Drummond and Hathrell [68] in a curved background 
context (see also Ref. [69]). Let us briefly review it in the CTP context. 

Functionally differentiating the CTP effective action we get the effective equations of 
motion for the expectation value of the field, ip := Tr(/)0) (see appendix [a| : 



6T[^pl,ip2] 



= 0. (5.4) 

Vl='fi2='fi 



6ipi{x 

The effective action can always be expanded in the proper vertices r"i ""''(xi, . . . ,Xr)- 

T[(fi,ip2]=^-^ I dV •••dVr"l-°'-(xi,...,X^r)<^ai(xi)---(^a,(Xr). (5.5) 

A straightforward generalization of the usual argument (see e.g. Ref. [49]) shows that 
this 2-point vertex corresponds to the inverse propagator, 

r'^\x,y) = i{G-y\x,y). (5.6) 

The Schwinger-Dyson equation, which defines the self-energy T,"'^{x,x'), 

Gab{x,y) = G^S{x,y)-i J d^z d^w G^^] {x, z)J:'\z,w)Gcb{w,y), 

can be manipulated to give 

{G-y\x, y) = A-\x, y) + iT.-\x, y), (5.7) 

where A"'^{x, y) is the inverse of the free propagator, 

A'^\x,y) = [(GW)-Y'(^,y) =c"M-rf + m2)5(4)(:c-y), (5.8) 

with 0°"^ = diag(l, —1). We see that the 2-point vertex can be expressed as 

T'^\x,y) = iA''\x,y) - T.-\x,y). (5.9) 

Hence the 2-point vertex essentially corresponds to the self-energy. Other proper vertices 
also have similar interpretations in terms of one-particle irreducible diagrams. 

If the relevant vertex is the 2-particle vertex T""^, the effective equations of motion can 
be expressed as 



5T 



d^y [T^\x, y) + T^\x, y)My) = 0, (5.10) 
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which, taking into account Eqs. (5.8) and (5.9) can be expanded as 



(_n^ + ^2^<^(^) + j d4ySR(x,y)99(y) =0, 



(5.11) 



where SR(x,y) = 'E^^{x,y) + E-^^(x,y) is the retarded self-energy. Introducing the 
Fourier transform around the mid point, the above equation can be rewritten as 



/ + ^2 + Sr(p; X) = 0. 



(5.12) 



Eq. (5.12) amounts to finding the poles of the retarded propagator. Thus, in flat space- 



time the self-energy and effective action methods lead to the same dispersion relation 
provided we use a CTP approach in both situations and we neglect vertices with three 
external particles or more. However, let us see that the interpretation of the dispersion 
relation is slightly different. Provided we can approximate. 



+ ^2 + Sr(p; X) 



-uj'^ +m'^ + Ep 



iujTr 



the solution to (|5.11|) can be written as 



(fix) 



d^k 



(2vr) = 



iEpt g-rpt/2 



(5.13) 



Therefore the real and imaginary parts of the poles have the role of the frequency of the 
damping rate of the mean field excitation. 



5.3. The mean-field and quasiparticle methods compared 

The mean-field-based approaches are simpler than the methods that we have developed 
in this paper, highlight the importance of the retarded propagator, and quickly relate 
the real and imaginary parts of the poles of the propagator to the energy and the decay 
rate respectively. However, the expectation value of the field operator is not a usual 
observable in field theory. Moreover, the kind of perturbations considered does not cor- 
respond to quasiparticles, since the latter, which are of the form a"^ pa, have vanishing 
expectation values for the field operator. Given all that, we believe that in relativistic 
field theory a linear response-based approach based on the study of the dynamics of ex- 
pectation value of the field operator is not well-suited to study the dynamical properties 
of quasiparticles in the regime of validity of the quasiparticle description. It is adequate 
though to describe the hydrodynamic regime. 

Relativistic field theory in the hydrodynamic limit [12, 13] can be understood as an 
effective theory describing the dynamics of long wavelength and short frequency per- 
turbations [14]. It has been extensively studied in the literature: for instance, let us 
mention that the evaluation of the hydrodynamic transport coefficients has been carried 
out by Jeon and Yaffe [70,71] and by Calzetta et al. [72], and that Aarts and Berges have 
studied the non-equilibrium time evolution of the spectral function in this regime [73] . 
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The limit of long wavelengths and short frequencies corresponds to considering small 
momenta. For Gaussian systems other than the vacuum, small-momentum modes have 
occupation numbers of the order of one or larger, an therefore this regime is not ad- 
equately described by a quasiparticle description, as commented in the introduction. 
Therefore, the hydrodynamic and quasiparticle descriptions represent two different com- 
plementary descriptions valid in two different regimes. The classical-like perturbations at 
low momenta are adequately described by a hydrodynamic description, and the individ- 
ual particle-like perturbations at high momenta can be analyzed within a quasiparticle 
formalism. 

In any case, we have seen that both the dynamics of the mean field and the quasi- 
particles can be characterized with the retarded self-energy. This can be understood 

in several ways, one of them being the following: loosely speaking, when the Gaussian 
approximation is introduced and higher order correlation functions are neglected, the 
system behaves as if it were effectively linear, and thus the dynamics of all relevant 
quantities is essentially determined by the solution of a corresponding stochastic prob- 
lem [21,45,74]. In this context it is not surprising that the elementary dynamics of both 
the mean field and the quasiparticles can be described with the same elements. 
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6. Summary and discussion 



In this paper we have presented two different approaches to the analysis of the quasipar- 
ticle properties in relativistic field theory: first, a frequency-based approach, in which 
the properties of the excitations have been deduced from the analysis of the spectral 
representation and, second, a real-time approach, wherein the quantum states corre- 
sponding to the quasiparticle excitations have been explicitly constructed, and the time 
evolution of the expectation value of the energy in those states has been studied. In both 
approaches it has been possible to show that the real and imaginary part of the self- 
energy determine the energy and the decay rate of the quasiparticles respectively — see 
Eqs. (3.8). Although previous evidences existed, to our knowledge this is the first sys- 
tematic corroboration that the real part of the self-energy determines the physical energy 
of quasiparticles in non-vacuum relativistic field theory. The dynamics of the quasipar- 
ticle excitation can be also encoded in the form of a generalized dispersion relation [see 
Eq. ^Hh^ y. 

£^ = El- iEpTp = m2 + p2 + Sr(£;p, p), (6.1) 

which needs not be Lorentz-invariant. 

Several additional interesting points have been illustrated by the real-time approach. 
We have explored the open quantum system viewpoint for the quantum field modes, 
and have found a physical criterion to fix the values of the system frequency and field 
renormalization parameters, equivalent to the on-shell renormalization condition in the 
vacuum. Using a very simplified model of a quasiparticle creation device, we have 
discussed that quasiparticles are adequately described by the action of the creation 
operator on the background state for long observation times — see Eq. (4.16). At the 
same time, we have shown that the quantum state corresponding to the quasiparticle 
contains actually more than one particle excitation. We have argued that this result is 
actually a statistical effect, due to the fact that sectors with higher occupation numbers 
are more likely to be excited, because of the Bose-Einstein statistics. We have also built 
the hole state [Eq. ( |4.9[ )], but have discussed that it cannot be properly considered a 
quasiparticle state for bosonic systems. 

Finally, by comparing our results with other mean- field-based approaches, we have 
argued that the latter are more suited for the study of classical-like configurations in the 
hydrodynamic limit rather than to the study of the properties of individual quasiparti- 
cles. 

Most expressions in this paper have referred to homogeneous, isotropic and stationary 
backgrounds. In fact, the main results depend, first, on the existence of the diago- 
nal relation (1.13) between the self-energy and the propagator in Fourier space (which 
demands homogeneity and stationarity); second, on the diagonalization of the density 
matrix in the basis of eigenstates of the Hamiltonian (which is implied by stationar- 
ity), and, third, on the equivalence of the p and — p modes (which requires isotropy). 
Therefore, it appears that the quasiparticle interpretation will be completely spoiled 
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when the background becomes non-stationary, non-homogeneous or non-isotropic. Let 
us argue that this is not the case, and that a particle interpretation is still feasible 
for non-homogeneous and non-stationary backgrounds provided that the characteristic 
scales of variation of the background state are sufficiently large. In the following we shall 
assume that L is the typical length scale in which the background changes significantly, 
and T is the typical timescale of evolution of the background. 

When the backgrounds are non-homogeneous or non-stationary there is still an ap- 
proximate diagonal relation between the retarded propagator and retarded self-energy 



[Eq. (1.13)], provided L ^ 1/E and T ^> 1/E, where E is the typical energy involved. 
The retarded propagator will have the usual analytic structure, with the only difference 
that now the location of the poles will be time and space dependent, and therefore the 
generalized dispersion relation will also be time and space dependent: 

£2=m^ + p^ + SR(^p(t, x), p; t, x) = 0. (6.2) 

On the other hand, for timescales much smaller than T, the background can be considered 
stationary. Hence, if one considers energies which correspond to those short timescales, 
E ^ l/r, the density matrix components corresponding to those energies are diagonal 
in the basis of eigenstates of the Hamiltonian. Therefore the properties which depend 
on the diagonalization, such as the spectral representation, are still valid. Finally, if the 
origin of the anisotropy is the inhomogeneity of the state, for the relevant energy scales 
the anisotropy will also be negligible if the condition L S> 1/E is fulfilled. However, it is 
possible that the system is homogeneous but anisotropic, and that the magnitude of the 
anisotropy is large. In this case some of the expressions we have written down will no 
longer be quantitatively valid, since we have demanded the equivalence of the forward 
and backward modes for a given momentum. However in Refs. [23, 26] it is shown with 
a particular example that the general picture is not essentially modified. 

Let us now comment the subtleties related to the construction of the quasiparticle 
states. Recall that quantum state corresponding to the quasiparticle contains more than 
one additional particle excitation. As we have already mentioned, this is a statistical 
effect related to the fact that the background state is not an eigenstate of the number 
operator. It can be illustrated with the following toy model. Let p = (|0)(0|-|-|2)(2|)/2 be 
the initial state of a harmonic oscillator. When a "quasiparticle" is introduced into the 
system by the action of the creation operator, the state becomes p' = (|1)(1| -|-3|3)(3|)/4. 
The expected number of particles in the initial and final states is 1 and 5/2 respectively, 
and therefore the particle number is increased in 3/2, despite the fact that each term 
of the state is increased with only one particle. Clearly, the reason for this is linked to 
the Bose-Einstein statistics: the higher occupation states are more likely to be excited. 
A similar phenomenon happens if the incoherent mixture is replaced by a coherent 
superposition of two particle eigenstates. See Ref. [75] for an analysis of the same 
phenomenon in another context, with a discussion of possible energy non-conservation 
issues. 

We have seen that the dispersion relations governing the particle dynamics can be 
also obtained from simpler mean-field-based methods (linear response theory or effective 
action), albeit with a slightly different interpretation. We believe that the effort put 
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in the construction and analysis of the quasiparticle excitations was anyway useful for 
different reasons. First, the quasiparticle and mean field approaches provide the same 
answer to two different physical questions: in the former case one studies the dynamics of 
individual quasiparticles, and in the latter, the dynamics of the mean fields. Second, and 
related to this, both descriptions have different complementary regimes of validity, the 
former approach being suited for high-momentum perturbations and the latter being 
suited for low-momentum perturbations. Finally, the route to the construction and 
analysis of the quasiparticle properties shed light on many intermediate results which 
are interesting by themselves. 

A key element in our analysis has been the Gaussian approximation, which on the one 
hand has allowed manageable expressions, and on the other hand has provided physi- 
cal interpretations for those expressions, without having to restore to any perturbative 
expansion in the coupling constant. The Gaussian truncation leads to the most basic 
description of the dynamics of the quasiparticles, the description in terms of a dispersion 
relation which we have investigated in this paper. Non-Gaussianities would appear in a 
more elaborate description of the quasiparticle dynamics. The Gaussian approximation 
can be formally controlled with the large- expansion in the number of fields [76,77]. 

In this paper we have limited ourselves to those systems which can be correctly de- 
scribed in terms of a pole in the retarded propagator, or, in other words, in terms a 
dispersion relation. As we have commented, there are systems posessing massless excita- 
tions in which the analytic structure of the propagator differs, leading to non-exponential 
decay laws [40-42,63] and to the so-called non-Fermi liquid behavior [64]. 

A natural extension of the work in this paper would be to generalize the results for 
non-zero spin fields. The analysis of fermion fields should reveal the emergence of hole 
excitations as true quasiparticles, in contrast to the bosonic systems, in which we have 
seen that hole states, although can be constructed, do not have the appropriate quasipar- 
ticle properties. Another direction in which the work of this paper could be extended is 
the analysis of quasiparticle interactions. In the same way as the Kallcn-Lehmann spec- 
tral representation can be adapted to non-vacuum situations, the Lehmann-Symanzik- 
Zimmerman formalism of vacuum field theory could in principle also be extended to 
non-vacuum situations. Some work has been already done in this direction in the con- 
text of thermofield dynamics [20] . 

The results of this paper can also be generalized to study the propagation of interacting 
particles or quasiparticles in curved backgrounds, the propagation in curved backgrounds 
having many similarities with the propagation in a physical medium [78,79]. In Ref. [65] 
the main results of this paper are generalized to account for the adiabatic propagation 
of interacting particles in cosmological backgrounds. 
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A. The closed time path method 



In this appendix we give a brief introduction to the closed time path (CTP) method 
(also called in-in method, in contrast to the conventional in-out method, or Keldysh- 
Schwinger method), stressing those aspects relevant for this paper. For further details 
we address the reader to Refs. [28,80-83]. 

Let us consider a free or an interacting scalar field (j). The path-ordered generating 
functional Zc[j] is defined as 

Zc[j] = Tt (pTc e^/cdt/d^x0(x)i(a;)\ ^ 



where (p^x) is the field operator in the Heisenberg picture, C is a certain path in the 
complex t plane, Tc means time ordering along this path and j(x) is a classical external 
source. By functional differentiation of the generating functional with respect to (p, 
path-ordered correlation functions can be obtained: 



(A.2) 

j=0 



Introducing a complete basis of eigenstates of the field operator in the Heisenberg picture, 
x)|(/), i) = (j){t,x)\(j),t) , as a representation of the identity, the generating functional 
can be expressed as: 

Zc[j] = [ d<Ad<A'(0,ti|/5|./.',ti)(0',ti|rce*/cd*/d^-'^(-)^-W|^,t.) (A.3) 



The functional measures d(f) and d(j)' go over all field configurations of the fields at fixed 
initial time t. If the path C begins and ends at the same point ti, then the transition 
element of the evolution operator can be computed via a path integral: 

Zc[j] = [d(j)d(j)'{^,t^\p\(l)',ti) r^*"''^ '^^''^p^gi/cdt/d3x{LM+^(x)iW}^ (^_4^ 

J ^93(ti,x)=(/>(x) 

where L[(j)\ is the Lagrangian density of the scalar field. 



Let us consider the time path shown in Fig. A.l, If we define (fi^2{t,^) = ip{t,x) and 



ji^2(i)X) = j(t,x) for t G Ci^2, then the generating functional can be reexpressed as: 
Z[ji,j2] = / d(t)d(t)'d(l)" {(p,ti\p\(l)' ,ti) 



X 



ipi(if,x)=(/)"{x) 

</3l{ti,x)=(/>(x) 



/•(^2(if,x)=(/)"(x) 
y</32{ti,x)=</.'{x) 
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Figure A.I.: Integration path in the complex-time plane used in the CTP method. The 
forward and backward lines are infinitesimally close to the real axis. 



In the following it will prove useful to use a condensed notation where neither the 
boundary conditions of the path integral nor the integrals over the initial and final 
states are explicit. With this simplified notation the above equation becomes 

Z[ji,j2] = J V(l)iV(P2{(p,t\p\(l)',t) e'f'^^''W'fii]-Hv2]+'fiiix)hix)-^2ix)j2ix)} (A_g) 
An operator representation is also possible: 

Z[ji,j2] = Tr T e"'4''^*-^'^^'''^('')-'2{3:) J. gi/t'/dt/d3x^(z)ji(x)^ ^ 

By functionally differentiating the generating functional the different correlation func- 
tions can be obtained: 



^ab\X,X ) - 



5j"-{x)5j^{x') 



ia=ib=o 



We recall that lowercase roman indices may acquire the values 1 and 2 are raised and 
lowered with the "CTP metric" Cab = diag(l,— 1). Higher order correlation functions 
can be obtained in a similar way. 

At least formally, all expressions so far are valid either for interacting or free field 
theories. However, explicit results can only be obtained when the theory is free and 



the initial state is Gaussian. In this case the path integrals in Eq. (A. 5) can be exactly 
performed, and one obtains: 

zW[ji,j2] = ZW[0,0] e-^/d^-d^-':''"(-)Gi^-,-')/(x') . (A.9) 

where the propagators G^^^ verify c°'^[—d^d'^ + m^jG^*^^ = —i6^'^\x — x'). As with the 
conventional in-out method, the perturbative expansion can be organized in terms of 
Feynman diagrams. There are two kinds of vertices, type 1 and type 2, and four kinds 
of propagators linking the two vertices. The Feynman rules are those of standard scalar 
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field theory, supplemented by the prescription of adding a minus sign for every type 2 
vertex. 

From the generating functional the connected generating functional j'] is defined 

W[j,j']:=-iliiZ[j,j'] (A.IO) 
Next we introduce the following objects: 

6W 

which must be understood as functionals of ji and j2 even if this dependence is not 
explicit. If ji = j2 both (pi and (p2 give the expectation value of the field under the 
presence of a classical source j. Finally, the effective action r[99i,(/32] is defined as the 
Legendre transform of the connected generating functional: 

r[v^i,^2] := W[h,j2] - I d^xja{x)^''{x). (A.12) 

In this equation ji and j2 must be understood as functionals of ipi and (p2, which can 



be obtained by inverting Eq. (A. 11). By functionally differentiating the effective action 
with respect to ifi and setting Lp2 = fi the equation of motion for the expectation value 
of the scalar field is obtained: 



6r 



6(p{x) 



= j{x). (A. 13) 

^2=^1=0 



In contrast to the conventional in-out treatment, the equations of motion obtained from 
the CTP generating functional are real and causal because they correspond to the dy- 
namics of true expectation values [84]. 

Thermal field theory can be seen as a strict particular case of the CTP method, in 
which the state p happens to be 



Tr(e 



-I3H] 



(A.14) 



where H is the Hamiltonian operator of the system. To apply the techniques presented in 
this appendix the only requirement is to compute the free thermal propagators. However, 
it proves convenient to make a slight a adaptation of the formalism and modify the 
complex time path in order to connect with the usual approaches to thermal field theory. 
Noticing that the density matrix operator can be seen as the time translation operator 
in the complex plane, e~^^ = U{t — if], t), Eq. (A. 3) can be reexpressed as: 



Zc[j]= / d0(</>,ti-i/3|Tce^/cd*/d^x0(.),(x)|^^^.^^ 

Tr(e~/^^) J 

where we have used the completeness relation J d(f)' \(f)' , t[) {(j)' , ti\ = 1. This way we have 
managed to incorporate the information on the state on the dynamical evolution. Now, 
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Figure A. 2.: Integration contour in the complex-time plane used in the real-time ap- 
proach to thermal field theory. The choice <t = 0^ makes the formalism 
analogous to the CTP method. 

if the path C starts at tj and ends at tj — if3, a path-integral representation can be 
introduced: 

Zc[j] = ^ / / p^gi/cdt/d3x{LM+^(z)i(x)}^ (^_^g^ 

Tr{e-f^H) J y</p(ti,x)=(/.(x) 

which can be rewritten in a more compact form as 

Zc[j]=N /"p(^e*i'cd*/d3x{LM+¥.(a;)i(x)}^ 



where is a normalization constant, and the boundary conditions 4>{ti, x) = (j){t\—i(3, x) 
are assumed. Different elections for the path C lead to different approaches to thermal 
field theory: a straight line from t\ to ti — 1(3 leads to the imaginary-time formalism, and 



the contour shown in Fig. A. 2 leads to the real-time formalism. By choosing (7 = 0^ the 
real-time formalism is virtually identical to the CTP method, since the the path along 
C3 and C4 can be neglected if we are interested in real-time correlation functions and the 
boundary conditions of the path integral are properly taken into account. 
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In the paper we make extensive use of the Gaussian states. In the following we shall 
give a brief description of the aspects relevant to us without entering into details. See 
e.g. Refs. [45,85,86] for a more complete introduction. In this appendix will always 
represent the correct normalization constant, which can always be computed with the 
Gaussian integration formula if desired. 

In general, Gaussian states are those whose density matrices have a Gaussian form 
in a coordinate representation, or equivalently those whose Wigner functions have a 
Gaussian form. A generic Gaussian state with (p) = and (q) = (or equivalently with 
(a) = (a^) = 0) can be represented as: 



P 



Nexp{- Fa^a + Gaa + G*a^a^), (B.l) 



where F is real. In this paper we have exclusively used Gaussian states with zero mean 
(many times without making explicit mention). Gaussian stationary states are those 
which conmute with the Hamiltonian, [p, H] = 0. The stationary state of a harmonic 
oscillator with H = $7 (a^^a + 1/2) must be of the form p = N exp {—Fa^a), and therefore 
corresponds to a thermal state. 

For a scalar field, decomposed in modes the most general state for the two-mode 
system itp which is translation invariant, i.e., which conmutes with the momentum 
operator]^ 

[/5s,p]=0, p = p(aj,ap - a;Lpa-p). (B.2) 
and which verifies (a±p) = (a±p) = is given by 

= iVexp [ - F(aJ,ap + alpfl-p) + 2Gapa_p + 2G*aJ,a;Lp] ■ (B.3) 
If we further impose stationarity with respect to the reduced Hamiltonian, 

[Ps: ^sys] = 0, -^sys = TTpTT-p + ^^pt^pC^-p = ^p(aj,ap + olpO-p + l) , (B.4) 

the most general Gaussian state corresponds to a factorized thermal state: 

Ps = iVexp [ - F(aJ,ap + a^pC-p)] = TVexp ( - FaJ,dp) exp ( - Fa!„pa„p). (B.5) 

A quantum mechanical or field theory system is Gaussian if its generating functional 
is Gaussian. For a single degree of freedom this means (assuming vanishing expectation 
values for the position operator) 



Z[ja] = exp 



(B.6) 



^Do not confuse the physical momentum operator p, with the canonical momentum operator TTp, 
conjugate of the field operator. 
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Gaussian systems correspond either to Gaussian states following quadratic equations of 
motion (as in Ref. [32]), or alternatively to an approximation for general states following 
general equations of motion (as in this paper). 

For Gaussian systems, according to the Wick theorem [6,8], the n-point correlation 
functions can be reduced to the two-point correlation functions. Instead of trying to 
give a general formulation of the Wick theorem let us simply show some particular 
applications. Whenever the system is Gaussian, any time-ordered four-point correlation 
function can be expressed in terms of two-point correlation functions (we assume that 
the expectation value of the field operators vanishes): 

(Tg(ti)g(t2)g(t3)g(t4)) = {Tqiti)q{t2)){Tq{t3)qit^)) + {Tq{ti)q{t3)){Tq{t2)q{t^)) 

+ {Tq{ti)q{ti)){Tq{t2)m)) 

If the correlation function is a mixture of time- and antitime-ordered expressions, the 
equivalent expression goes as follows 

{Tq{ti)q{t2)fq{t3)q{ti)) = {Tq{ti)q{t2)){fq{t3)q{t^)) + (g(ti)g(t3)) (^(^2)9(^3)) 

+ mi)qit,)){q{t2)m) 

These expressions can be demonstrated by taking derivatives on the Gaussian generating 
functional (B.6). 

In field theory, if the state is translation-invariant, momentum conservation can sim- 
plify the application of the Wick theorem. For instance, if k 7^ q, 

{T4>k{tl)4>-ci{t2)Mt3H-UU)) = {TMh)4>-Ut4)) + (r0q(t2)0-q(t3)). (B.9) 

The two-point correlators can be also expressed as a function of the creation and anni- 
hilation operator. If the state is stationary, only those terms having a creation and an 
annihilation operator survive: 

{Mti)4>-Ui2)) = ^ ((a^_k(ti)a_k(i2)) + (ak(ii)4(*2))) • (B.IO) 
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C. Contour integration of I{t,to;p) 



When computing the time evolution of the propagator, the fohowing integral appears 



see eqs. (4.23b) and (4.23c)]: 



Apparently, this integral cannot be evaluated with complex plane techniques since the 
integrand contains the factor and |a;|, which is non-analytic with the usual prescription 
= \/iOio*, and the factor 9 {to), which in principle is defined in the real axis. Let us 
do a more careful analysis. 

We begin by extending the problematic terms to the complex plane in the following 
way: 

9{ll!) 6{Reu)), \uj\ w sign(ReLij). 

Those terms continue to be non-analytic, but only on a branch cut located in the imag- 
inary axis. Therefore, with this prescription the integrand is analytic everywhere on the 
complex plane except on the branch cut and on the poles. Notice that the branch cut is 
rather special, since the function is continuous across the branch cut. 

Second, notice that the following related contour integrals are well-defined for t > Iq 

hit, to; P) = — / - e--(*-*o) -^(^ + ^p)^rpnp 

Uit, to; p) = 4. / ^ e--(*-o) + ^p)'rpK + 1) 



G4 

where C3 and C4 are the closed anticlockwise paths at the boundaries of the left lower 
and right lower quadrants respectively 

Third, note that with the above prescription /(t, to; p) = /3(t, to; p) + hit, to; p), since 
the path at infinity does not contribute, and the contribution from the path in the 
imaginary axis cancels (because the function is continuous at the branch singularity). 
Therefore, in practice, the integral /(t, to;p) can be computed by residues as if there 
were no branch cut singularity. 

Let us now compute the integrals l3(t, to;p) and /4(t, to;p). We start by l4(t, to;p). 
There is a pole in the lower right quadrant at a; ~ Ep — iT p/2. Neglecting Fp in front 
of Ep we find 

l4(t, to; p) « Epil + Up) e-^^p(*-*o) e-rp/2 _ 

With respect to the integral /4(t,to;p), the contribution from the pole in the lower left 
quadrant is very suppressed because of the factor iuj-\-Ep)'^ in the numerator. Therefore 
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/3(t, to; p) ^ -^4(^5 ^o; p)- Given all this we obtain the final result: 

I{t, to; p) ~ Ep{l + np) e-^-^p(*-*o) e-^p/^ . (C.2) 
Alternatively this integral can also be evaluated directly in the time domain [65] . 
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